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Gravitational waves from cosmic strings are generated in the first fractions of a 
second after the Big Bang, potentially providing a unprecedented probe of the early 
universe. We discuss the key dynamical processes underlying calculations of the 
stochastic background produced by a string network and we detail the parameter 
dependencies of the resulting spectral density fl gT (f). The present constraints 
on the cosmic string mass scale // arising from the millisecond pulsar timings and 
primordial nucleosynthesis are discussed, justifying our present conservative bound 
Gp/c 2 < 5.4(±1.1) X 10 C We then discuss the strong prospect of detecting (or 
ruling out) cosmic string background with the next generation of gravitational 
wave experiments. Comparison is also made with alternative cosmological sources 
of gravitational waves such as inflation and hybrid topological defects. 


1 Introduction 

An outstanding achievement of modern cosmology, has been the detection of 
anisotropies in the cosmic microwave background!!]. These anisotropies pro¬ 
vide a snapshot of the universe about 400,000 years after the Hot Big Bang, 
just as the universe became transparent to electromagnetic radiation. Despite 
this observational triumph, as yet it has been insufficient to differentiate be¬ 
tween competing paradigms for galaxy formation, nor has it shed much light 
on cosmological processes taking place before photon decoupling. There are, 
however, other types of radiation. Gravitational radiation, as considered in 
the present work, may penetrate through this electromagnetic surface of last 
scattering, and travel virtually unaffected since emission. Of course, for gravi¬ 
tons this remarkable transparency is due to their very weak interaction with 
ordinary matter which, in turn, makes them difficult to observe. However, 
pioneering experiments have been proposed which could detect a stochastic 
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background of gravitational waves generated in the early universe over a range 
of frequencies. 

Cosmic strings are line-like topological defects which may have formed 
during a phase transition in the early universeota. Strings which formed with 
a mass-per-unit-length /x such that G/x/c 2 ~ 10~ 6 may be responsible for 
the formation of the large-scale structure and cosmic microwave background 
anisotropy observed in the universe today. In general, a network of strings 
will evolve towards a self-similar scaling regime by the production of loops 
and subsequent emission of gravitational radiation. Since the sources of these 
gravitational waves come from many horizon volumes this background will be 
stochastic and is likely to appear as noise on a gravitational wave detector. 

A stochastic background of gravitational waves is normally quantified by 
the relative spectral density fl g (/) given at a frequency /. That is, the en¬ 
ergy density in gravitational radiation in an octave frequency bin centred on 
/, relative to the critical density of the universe. This is directly related to the 
dimensionless wave amplitude ( h c oc f) which is measured experimen¬ 

tally. 

In the case of gravitational radiation produced by a cosmic string network, 
for a given frequency today we may identify a characteristic time at which the 
waves were emitted. Assuming that the radiation is emitted at a time t e before 
equal matter-radiation (t e < t eq ~ 4, 000 years) and that it is created with a 
wavelength comparable to the horizon A (t e ) ~ t e , then the frequency today 
is given by / ~ Aiq^eq^e) -1 / 2 where the redshift is z eq ~ 2.3 x 10 4 rio^ 2 - In 
particular those gravitational waves created at the electroweak phase transi¬ 
tion will have a frequency / ~ 10 _3 Hz, whereas those created at the time of 
radiation-matter equality will have a frequency / ~ 10 _13 Hz. 

The relevant constraints on stochastic backgrounds are due to the mea¬ 
sured timing residuals in pulsar signals and also from Big Bang nucleosynthesis. 
The most recent analysis of pulsar signal arrival times gives the limit on the 
spectral density of gravitational radiation 

ngr(/) = I /oh. < 9-3 X 10- 8 *.- 2 (95% CL) (1) 

Per it 

in a logarithmic frequency interval at / 0 b s = (8yrs) -1 . This analysis corrects 
errors in previous workHQ, and uses an improved method for testing the hypoth¬ 
esis that the timing noise is due to a stochastic gravitational wave background. 
For these reasons we use the latest bound, equation (Q), to obtain a constraint 
on G/x/c 2 for given values of the cosmic string and cosmological parameters. 

The bound from nucleosynthesis is slightly different, as it constrains the 
entire spectrum. The standard constraint quoted is that the total energy den- 
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sity in gravitational waves at the time of nucleosynthesis must be less than 
5.4% of the critical density, or 

/ 7 

— f l gr (u;) « 0.162 (N v - 3)fl r < 0.054, (2) 

LO 

which can be thought of as a contribution from an extra neutrino species. 
This constraint is somewhat weakened at present, due to uncertainties in the 
observed element abundances. For this reason, we shall discuss constraints 
from a range of values for N„, the effective number of neutrino species. 

In this proceedings, based in part on work in ref. [7], we review the status 
of the stochastic background produced by cosmic strings. The details of the 
spectrum of gravitational radiation due to cosmic strings has been previously 
considered elsewhereO. Here we benefit from recent works, which suggests that 
the effect of the radiative back-reaction on a string loop is to damp out the 
higher oscillation modes. We illustrate the effects of a non-standard thermal 
history, low-fl universes and discuss hybrid systems of defects, such as strings 
and domain walls or strings connected by monopoles. Finally, we discuss the 
opportunities for detection. 


2 Cosmic string loop radiation spectrum 

The spectrum of gravitational radiation emitted by a network of cosmic string 
loops is obtained using a background Q = 1 FRW cosmological model, an 
extended one-scale model for the evolution of a network of cosmic strings, and 
a model of the emission of gravitational radiation by cosmic string loops. The 
procedure by which.the spectrum is computed has been presented in detail by 
Caldwell and AllenEI. In this section we discuss the model for radiation by an 
individual loop. 

The model of the emission of gravitational radiation by cosmic string loops 
is composed of the following three elements. 

1. A loop radiates with power P = TG^ 2 c. The dimensionless radiation 
efficiency, T, depends only on the loop-configuration, rather than overall 
size. Recent studies of realistic loopsll3 indicate that the distribution of 
values of the efficiency has a mean value (F) ss 60. 

2. The frequency of radiation emitted by a loop of invariant length L is 
f n = 2 n/L where n = 1,2,3,... labels the oscillation mode. 
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3. The fraction of the total power emitted in each mode of oscillation n at 
frequency f n is given by the coefficient P n where 

OO 

p=(J2 Pn ) G ^ c = TG ^ c - ( 3 ) 

n =1 

Analytic and numerical studies suggest that the radiation efficiency co¬ 
efficients behave as P n oc n~ q where q is the spectral index. 

This model has several shortcomings. First, the spectral index q has not been 
well determined by the numerical simulations. Numerical workO suggests 
q = 4/3 as occurs with cuspy loops, loops along which points momentarily 
reach the velocity of light, based on simulations of a network of cosmic strings. 
However, these simulations have limited resolution of the important small scale 
features of the long stringpand loops. Hence, the evidence for q = 4/3 is not 
compelling. Analytic worklla suggests that q = 2, characteristic of kinky loops, 
loops along which the tangent vector changes discontinuously as a result of in- 
tercomnrutation, may be more realistic. Second, the effect of back-reaction on 
the motion of the string has been ignored. In this model, a loop radiates at all 
times with a fixed efficiency, in all modes, until the loop vanishes. As we shall 
next argue, the back-reaction will result in an effective high frequency cut-off 
in the oscillation mode number. Thus, the loop will only radiate in a finite 
number of modes, and hence in a finite range of frequencies. The resolution of 
these issues may have strong consequences for the entire spectrum produced 
by a network of strings. 

Recent advances in the understanding of radiation back-reaction on global 
strings suggest various modifications to the simplified model of emission by 
cosmic string loops. It has been shownEj that there are remarkable similarities 
between gravitational radiation and Goldstone boson radiation from strings, 
which we believe allow us to make strong inferences as to the nature of grav¬ 
itational radiation back-reaction. For example, the same, simple model for 
the emission of gravitational radiation by cosmic strings may be transferred 
over to global strings: in the absence of Goldstone back-reaction, global string 
loops radiate at a constant rate, at wavelengths given by even sub-multiples 
of the loop length, with an efficiency as described by an equation similar to 
(§). Hence, our argument proceeds as follows. Fully relativistic field theory 
simulations of global strings have been carried outd, where it was observed that 
the power in high oscillation modes is damped by the Goldstone back-reaction 
on periodic global strings. An analytic model of Goldstone back-reactionEl, as 
a modification of the classical Nambu-Goto equations of motion for string, was 
developed which successfully reproduces the behaviour observed in field theory 
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simulations. That is, high frequency modes are damped rapidly, whereas low 
frequency modes are not. Thus, we are motivated to rewrite equation (|^) for 
global strings, and by analogy for cosmic strings, as 


P= (j2Pn)G^c = TG^ 2 


( 4 ) 


n=1 


where n* is a cut-offfi-3 introduced to incorporate the effects of back-reaction. 
By comparing the back-reaction length-scale to the loop size, we estimate 
that such a cut-off should be no larger than ~ (TGuVc 2 ) -1 . The ongoing 
investigations of global and cosmic string back-reactiontj have not yet reached 
the level of precision where a firm value of n* may be given. As we demonstrate 
later, the effect on the radiation spectrum is significant only for certain values 
of the cut-off. 

3 Analytic estimate of the radiation spectrum 

Analytic expressions for the spectrum of gravitational xadiation emitted by a 
network of cosmic strings have been derived elsewhereB’0. While these analytic 
expressions are simplified for convenience, they offer the opportunity to exam¬ 
ine the various dependencies of the spectrum on cosmic string and cosmological 
parameters. 

The spectrum of gravitational radiation produced by a network of cosmic 
strings has two main features. First is the ‘red noise’ portion of the spec¬ 
trum with nearly equal gravitational radiation energy density per logarithmic 
frequency interval, spanning the frequency range 10“ 8 Hzl$/;$10 lo Hz. This 
spectrum corresponds to gravitational waves emitted during the radiation- 
dominated expansion era. This feature of the spectrum may be accessible to 
the forthcoming generation of gravitational wave detectors. Second is the peak 
in the spectrum near / ~ 10“ 12 Hz. The amplitude and slope of the spectrum 
from the peak down to the flat portion of the spectrum is tightly constrained 
by the observed limits on pulsar timing noise. 

3.1 Red noise portion of the spectrum 

An analytic expression for the ‘red noise’ portion of the gravitational wave 
spectrum is given as follows: 




( 5 ) 
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A = Pood 2 H (t)c 2 / g (3 = [1 + f T ad H {t)c/(TGgt)\ \ (6) 

In the above expressions, poo is the energy density in ‘infinite’ or long cosmic 
strings, a is the invariant length of a loop as a fraction of the physical horizon 
length dn{t) at the time of formation, (v 2 ) is the rms velocity of the long 
strings, and / r « 0.7 is a correction for the damping of the relativistic center- 
of-mass velocity of newly formed string loops. All quantities are evaluated in 
the radiation era; dn{t) = 2ct, A ==.252 ± 10, and ( v 2 )/c 2 = 0.43 ± 0.02 as 
obtained from numerical simulational3. 

The above expression for the spectrum has been obtained assuming no 
change in the number of relativistic degrees of freedom, g, of the background 
radiation-dominated fluid. However, the annihilation of massive particle species 
as the cosmological fluid cools leads to a decrease in the number of degrees of 
freedom, and a redshifting of all relativistic particles not thermally coupled 
to the jflaiid. This has the effect of modifying the amplitude of the spectral 
densityEJ, equation (|j), by a factor (g(Tf) /g(Tj)) 1 / 3 where g(Tij) is the num¬ 
ber of degrees of freedom at temperatures before and after the annihilations. 
Using a minimal GUT particle physics model as the basis of the standard ther¬ 
mal history, we see that the red noise spectrum steps downwards with growing 
frequency. 


f dp gr 

Pcrit df 


-1) 

9 ac 4 (z e q +1) 


1 


10 -8 Hz < / < 10 -10 O! -1 Hz 


(3.36/10.75) 1 / 3 = 0.68 
X< 10 -10 a -1 Hz < / < 10 -4 a -1 Hz 

(3.36/106.75) 1 / 3 = 0.32 

10 -4 a -1 Hz < / < 10 8 Hz 


(7) 


Hence, the red noise spectrum is sensitive to the thermal history of the cosmo¬ 
logical fluid. The locations of the steps in the spectrum are determined by the 
number of relativistic degrees of freedom as a function of temperature, g{T). 
As an example, we present the effect of a non-standard thermal history on the 
spectrum in Figure In this sample model, the number of degrees of freedom 
g{T) decreases by a factor of 10 at the temperatures T = 10 9 , 10 5 , 1 GeV. The 
effect on the spectrum is a series of steps down in amplitude with increasing 
frequency; detection of such a shift would provide unique insight into the par¬ 
ticle physics content of the early universe at temperatures much higher than 
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Figure 1: The effect of a non-standard thermal history of the cosmological fluid on the 
amplitude of the red noise portion of the gravitational wave spectrum is shown. The solid 
curve displays the spectrum produced using a minimal GUT with a maximum g = 106.75. 
The dashed curve shows the spectrum produced allowing for a hypothetical, non-standard 
evolution of g(T), as might occur if there were a series of phase transitions, or a number 
of massive particle annihilations as the universe cooled. For temperatures T > 10 9 GeV, 
the number of degrees of freedom is g = 10 4 . For 10 5 GeV < T < 10 9 GeV, g = 10 3 . For 
T < 10 5 GeV, the standard thermal scenario is resumed. 


may be achieved by terrestrial particle accelerators. In the case of a cosmo¬ 
logical model with a thermal history such that g(T) g(Tf) for X > Tf , all 
radiation emitted before the cosmological fluid cools to X will be redshifted 
away by the time the fluid reaches Tf. As we discuss later, such a sensitivity 
of the spectrum to the thermal history affects the nucleosynthesis bound on 
the total energy in gravitational radiation, and the opportunity to detect high 
frequency gravitational waves. 

3.2 Peaked portion of the spectrum 

We now turn our attention to the peaked portion of the gravitational wave 
spectrum. The shape of this portion depends critically on the model for the 
emission by a loop, presented in section |^. Specifically, if a loop emits in very 
high modes, a significant portion of the loop energy maybe radiated at high 
frequencies, well above the fundamental frequency of the loop. Hence, the 
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Figure 2: The effect of a cut-off in the radiation mode number on the spectrum of gravita¬ 
tional radiation is shown. Curves for the loop radiation spectral index q = 2, 4/3 for various 
values of n* are shown. The vertical line shows the location of the frequency bin probed by 
pulsar timing measurements. For n* < 10 2 the shape of the spectrum is insensitive to the 
value of q for purposes of pulsar timing measurements. For increasing n*, more radiation 
due to late-time cosmic string loops is emitted in the pulsar timing frequency band. 


dominant behaviour of the peaked portion of the spectrum is given by 

J_ dpgr_ _ f Ci//^- 1 ) 1 < q < 2 

Pcrit df ~ \ C 2 /f q> 2 

for 10" 12 Hz</<10" 8 Hz. (8) 

Here GY 2 are dimensionful quantities which depend on G/r/c, a, T, A, q and 
n*. A lengthy expression displaying the full dependence of the spectrum 
on these parameters is not particularly enlightening, they can be found in 
ref.[13].The qualitative behaviour is as follows. The overall height of the spec¬ 
trum depends linearly on G/r/c 2 , while the frequency at which the peaked 
spectrum gives way to the red noise spectrum depends inversely on a. The 
important result is that for values of the mode cut-off n*<10 2 , the spectrum 
drops off as 1 // for any value q > 4/3. As a demonstration, sample spectra 
with various values of n* are displayed in Figure || Hence, the introduction of 
a sufficiently low mode cut-off eliminates the dependence of the spectrum on 
q , the loop spectral index. 








3.3 How stochastic is the background? 

We now comment on the statistical properties of the background. We would 
like to know how good an approximation it is to claim that the signal is stochas¬ 
tic. The relevant quantity to consider is the number of string loops which con¬ 
tribute to the radiation produced in a particular frequency bin. If we identify, 
for each frequency / observed, a unique time of emission corresponding to the 
time tf that the loop formed, then one can calculate the number of individ¬ 
ual horizon cells on the sky contributing to the radiation background at each 
frequency: 

N(f) « 4 x 10 1 'V (£) . (9) 

Here, f p « 4 x 10 _9 Hz is the frequency corresponding to horizon sized grav¬ 
itational waves that would be detected by the milli-second pulsar and a is 
the loop production size with respect to the horizon. For the frequencies of 
interest, / « 100Hz (LIGO/VIRGO), / « 10" 3 Hz (LISA) and / w 4 x 10" 9 Hz 
(pulsar), the values of N are 10 28 , 10 18 and 4 x 10' respectively, where we have 
used a ~ 10 -5 for illustrative purposes. Hence, we may argue from the central 
limit theorem, that if the N cells act independently, the sum of the amplitudes 
of radiation arriving at the detector or pulsar is gaussian and the background 
is stochastic. This gaussian assumption may break down at a frequency of 
10~ 12 Hz where TV ~ 1, corresponding to radiation emitted around a redshift 
z ~ 4. In this case, the radiation is due to a relatively small number of sources 
which are on average a few hundred megaparsecs away. These sources would 
appear as burst sources at a detector, rather than continuous noise. Unfor¬ 
tunately, this region of the spectrum is outside the range of interferometers 
since the timescales (~ 10,000 years) required for a detection are too large, 
but it may be possible to detect these gravitational waves using, for example, 
gravitational lensing. 

4 Alternative scenarios 

4-1 Open universes 

There is a substantial body of astronomical evidence which suggests that the 
cosmological density parameter is less than critical, il < 1. Therefore, it seems 
sensible to consider the effect of a low-density universe on the spectrum of 
gravitational radiation. 

The evolution of strings in an open universe will be very much the same 
as for the flat case, except that after curvature domination at t curv ~ Oot-o, 
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Figure 3: The effect of a low density, Qo < 1 universe on the peaked portion of the 
gravitational wave spectrum. The solid, long- and short-dashed curves represent spectra for 
Slo = 1, 0.6, 0.2. The vertical line shows the location of the frequency bin probed by pulsar 
timing measurements. For the loop spectral index q = 2, a low density universe dilutes only 
the lowest frequency waves, corresponding the radiation emitted by loops still present today. 


where a(t) jx the linear scaling regime no longer exists. In fact, it has 
been shownta that the characteristic length-scale of the network will increase 
like L ~ t[log(t)] 1//2 , rather than L ~ t, and hence the long string density 
Poo oc L~ 2 oc [t 2 log(f)] -1 , decreases relative the critical density p cr it oc t~ 2 , 
but increases with respect to the background density po oc t~ 3 . This is impor¬ 
tant for consideration of the gravitational waves created with low frequencies 
in the matter era and also for the normalization of the cosmic microwave back¬ 
grounds] as discussed in section ||. However, from the point of view of the 
red noise spectrum, which is produced in the radiation dominated era, there 
is little difference apart from a shift in the frequency corresponding to equal 
matter-radiation. 

Using methods similar to those used in the previous section, but modified 
to accommodate < 1, we have examined the spectrum of gravitational radi¬ 
ation produced by the cosmic string network in an open FRW space-time, with 
0.1 < f2o < 1. For this range of values of the cosmological density parameter, 
the portion of the spectrum produced at a time t is shifted downward by a 
factor ~ f }(t). Sample spectra for various values of Qq are displayed in Figures 
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Figure 4: The effect of a low density, Qo < 1 universe on the peaked portion of the 
gravitational wave spectrum. The solid, long- and short-dashed curves represent spectra for 
f2o = 1, 0.6, 0.2. The vertical line shows the location of the frequency bin probed by pulsar 
timing measurements. For the loop spectral index q = 4/3, a low density universe leads to 
a dilution of gravitational waves with wavelengths up to / 10 —5 Hz. 


||||. In the case that the spectrum drops off slower than 1//, that is 1 < q < 2 
and n* —► oo, the spectral density at frequencies as high as / ~ 10~ 5 Hz is 
diluted for Q 0 < 1. In the case that the spectrum drops off as 1//, that is 
q > 2, only at lower frequencies, /<10 _lo Hz, is the spectral density affected. 
This behaviour is exactly as expected: only those frequencies dominated by 
matter eras loops are affected. 

4-2 Hybrid defects 

Cosmic strings are generic in a number of symmetry breaking schemes where 
the first homotopy group is non-trivial. Simple examples where this can take 
place are, when the original symmetry group contains a U( 1) subgroup which is 
broken to the identity or when the broken symmetry group has a Zi subgroup. 
However, if these simple examples occur as part of a much more complicated 
scheme, it is likely that hybrid systems such as strings connected to either 
domain walls or monopoles may form. In the case when the hybrid system 
does not annihilate immediately this may lead to a stochastic background in 
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a similar way to that for stringar9. We will, therefore, illustrate the basic 
physics of this possibility by reference to the two examples already discussed. 
The interesting feature of these models being that they evade the constraints 
from the cosmic microwave background and pulsar timing allowing for larger 
values of G/z/c 2 and hence larger contributions to the stochastic background 
in the detectable range of frequencies. 

The first and simplest case is that of domain walls connected to strings. 
This will occur in the following symmetry breaking scheme, 

G-^H x Z 2 -+H, (10) 

where at the first transition (T = r/ s ) strings form and at the second (T = ?y w ) 
each string gets attached to a domain wall. Assuming that the domain walls 
form during the radiation era, then before the formation of walls the strings 
will evolve as in the standard scenario, creating a nearly flat stochastic back¬ 
ground in the frequency range / £ [2/a(t w t eq ) 1 ^ 2 , 2/a(f*f e q) 1 ^ 2 ], where td = 
min [fx/a,t w \ is the time at which the walls dominate the dynamics of the net¬ 
work, f w is the time of wall formation, f * is the time when relativistic evolution 
of the string network begins and /i ~ 77J, cr ~ rf^ are the string tension and wall 
surface tension respectively. Once the walls dominate the dynamics of the net¬ 
work, it will break up into what are effectively string loops spanned by domain 
walls. These string loops will collapse into gravitational radiation and other 
decay products, such as gauge particles, in about a Hubble time. The exact 
nature of this contribution is obviously dependent on the distribution of the 
loops produced by the fragmentation process. This process takes place over a 
relatively short frequency range and may lead to a sharp spike in the spectrum. 
If a sharp peak in the background spectrum was detected, the the demise of 
hybrid defects might provide one potential explanation. However, unlike the 
flat string background, the sensitive dependence on model phenomenology in 
this case, means that it is not possible at present to predict whether or where 
such a peak should occur. 

The second possibility is of strings connected by monopoles which would 
have been created in a symmetry breaking scheme such as, 

G -> H x 17(1) -> H , (11) 

where monopoles are formed at the first transition (T = ?7 m ) and strings 
connected between monopole/anti-monopole pairs are formed at the second 
(T = 77s). In this case there is no period in which the strings behave as they do 
in the standard scenario. If the first transition is after any period of inflation 
which may have taken place, then the average separation of the monopoles 
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is less the Hubble radius and the hybrid system will collapse in one Hubble 
time by dissipating energy into friction with the cosmological fluid. However, 
if the monopoles are formed during inflation then the system does not collapse 
immediately, allowing theibrmation of a stochastic background. This scenario 
has been studied in detaitH for the simplest case of a straight string connect¬ 
ing two monopoles. It was found that the power spectrum emitted by this 
simple configuration is divergent with P n oc n -1 due to the perfect symmetry 
(a similar spectrum is emitted by a perfectly circular string loop). Assuming 
that this divergence is softened in a less symmetric case, one may be able to 
calculate accurately the stochastic background created. However, it is likely 
to be dependent on the initial distribution of monopoles created in the first 
phase transition. 


5 Observational bounds on the radiation spectrum 


In this section we determine the observational constraint on the cosmic string 
mass-per-unit-length G/i/c 2 . To begin, we discuss the recent analyses of the 
pulsar timing data, after which we apply the newly obtained bounds to the 
cosmic string gravitational wave background. 

The observations used to place a limit on the amplitude of a stochas¬ 
tic gravitational wave background^ consist of pulse arrival times for PSR 
B1937+21 and PSR B1855+09. Although there has been some recent contro¬ 
versy^ regarding the analysis of the pulsar timing dataQ, the work by McHugh 
et acl best assesses the likelihood that the timing residuals are due to grav¬ 
itational radiation. We note that all analyses to date have assumed a flat, 
red noise spectrum for the gravitational wave spectral density. Such an as¬ 
sumption is only justified for a restricted range of frequencies in the case of a 
background due to cosmic strings, as we have demonstrated in the preceding 
section. Hence, a statistical analysis which uses a realistic model of the cosmic 
string spectrum may obtain a different limit on the amplitude of the spectral 
density. 

We now present values of the parameter G/x/c 2 for values of a which satisfy 
the pulsar timing constraint on the gravitational radiation spectrum. Contours 
of constant fl gr in the logarithmic frequency bin / = (8yrs) —1 , given by , in 
(a, G/Lt/c 2 ) parameter space, are shown in Figure [|. We have used cosmological 
parameters ff 0 = 1 and h G [0.5, 0.75] with the cosmic string loop radiation 
efficiency T = 60. We find 


. , f 2.0(±0.4) x 10- 6 (2/i) -8 / 3 q = 4/3 

\5.4(±1.1) x 10“ 6 q > 2 or n*<10 2 . 


( 12 ) 
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Figure 5: Curves of constant f2 gr in (a, Cji/c 2 ) parameter space are shown. For a given 
value of a, these figures give the observational bound on Gp/c 2 in the case h = 0.5, 0.75. 
In each figure, the constraining curves for g = 10, 2, 4/3 are given by the solid, long-, and 
short-dashed curves. The light dashed lines show a = TG/r/c 2 . The most conservative 
constraint is Gp/c 2 < 5.4 X 10 -6 . 


These constraints correspond to the maximum value of Gfi/c 2 along the con¬ 
tour of constant fl gr . In the case q = 4/3, this maximum occurs near a = 
TG/d/c 2 , the expected size of newly formed loops based on considerations of 
the gravitational back-reaction, while for the q > 2 or n* < 10 2 case, the maxi¬ 
mum occurs at a slightly smaller value of a. For both larger and smaller values 
of a the bounds become more stringentn. The unusual dependence on h is due 
to the contribution from high mode number waves emitted in the matter era, 
for which the amplitude depends on both the slope of the spectrum and the 
time of radiation-matter equality. The quoted errors are due to uncertainties_m 
the cosmic string model parameters measured by the numerical simulational3. 
For the case of an open universe, there is no change in the q > 2 or n* < 10 2 
bound. However, the q = 4/3 bound is weakened by a factor ~ 1/fio- 

We now comment on the validity of the model which we have used to 
generate the gravitational radiation spectra. We have shown that the obser¬ 
vational bounds on the total spectrum are sensitive to the value of the loop 
spectral index q , unless there is a back-reaction cut-off n* ^ 10 2 . Furthermore, 
we have noted that there is uncertainty in the characteristic value of the loop 
spectral index, q. Hence, we feel that it is more reasonable to take the con¬ 
servative bound of (^||) at the present. Next, consider the extended one-scale 
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modeKra, for the evolution of the string network. This model assumes that 
the long string energy density scales relative to the background energy density, 
with the dominant energy loss mechanism due to the formation o£ loops of a 
characteristic scale. A more sophisticated model, by Austin et al E3, attempts 
to include the effect of the gravitational back-reaction on the long-term evolu¬ 
tion of the string network; results suggest that an effect of the back-reaction 
may be to lower the scaling density in long strings at late times, beyond the 
reach of numerical simulations. Hence, there is some uncertainty as to how 
accurately the extended one-scale model describes the evolution of the string 
network. However, we do not believe that these considerations could result in 
a decrease in the amplitude of the gravitational wave background by more than 
~ 50%. Thus, we quote G/x/c 2 < 5.4(±1.1) x 10~ 6 as a conservative bound on 
the cosmic string mass-per-unit-length. 

The bounds already computed!! due to the constraint on the total energy 
density in gravitational waves at the time of nucleosynthesis remain valid. For 
a limit on the effective number of neutrino species N v < 3.1, 3.3, 3.6, the 
bound on the cosmic string mass-per-unit-length is G/x/c 2 < 2, 6, 10 x 1CU 6 
respectively, evaluated at a = FG/x/c 2 . The big-bang nucleosynthesis limit 
on the number of effective neutrino species is a conservative 7V„ < 4, owing 
to uncertainties in the systematic errors in the observations of light element 
abundancescn. Hence, until the observations are refined, the nucleosynthesis 
bound is weaker than the pulsar timing bound. Furthermore, the translation 
of the limit on N v into the_bound on G/x/c 2 is sensitive to the thermal history 
of the cosmological fluicfiZl. The bound on the string mass-per-unit-length 
may be considerably weakened if the cosmological fluid possessed many more 
relativistic degrees of freedom in the early universe beyond those given by a 
minimal GUT model. 

Comparing detailed computations of the large angular scale cosmic mi¬ 
crowave background temperature anisotropies induced by cosmic stringaij with 
observations, the cosmic string mass-per-unit-length has been normalized to 


G^/c 2 = 1.05" x 1CT 6 . (13) 

In an open universe, this normalization is expected to increase slightlyEl to 
G/x/c 2 ~ 1.6 x 10 -6 for = 0.1. Therefore, given the uncertainties in the 
extended one-scale model, we find the gravitational radiation spectrum to be 
compatible with observations. 

We have already noted that the hybrid systems evade the constraint from 
pulsar timing unless they survive until around the time of nucleosynthesis 
(tnuci ~ lsec), much later than thought possible according to the standard 
model of particle physics. Therefore, the only constraint on the stochastic 
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background comes from Big Bang nucleosynthesis. The interesting thing to 
note is that since the hybrid network will annihilate at some time f a , well 
before nucleosynthesis, then assuming that the spectrum is flat, the amplitude 
of the stochastic background allowed is larger by a factor of 

log (¥)/ log G:)- <14) 

For the first transition at the GUT scale, that is f* ~ 10~ 32 , and the second 
around 10 lo GeV with almost immediate annihilation, then this results in the 
constraint being weakened by approximately a factor of 10, corresponding to 
G/x ~ 10 -5 . This makes a detection in the first generation of LIGO detectors 
marginally possiblec2l, although the parameter values would have to be rather 
extreme. 

6 Detection of the Radiation Spectrum 

We would like to determine whether the stochastic gravitational wave spectrum 
emitted by cosmic strings may be observed by current and planned detectors. 
Because all ground-based detectors operate at frequencies /> 10 -3 Hz, we need 
only consider the ‘red noise’ portion of the gravitational wave spectrum (71). 
Noting that the spectral density, Sl gr (/), has a minimum value when a —> 0 
the predicted spectrum is bounded from below by 

*W) > wrA/r (1 ~ j^ } c2) (Gp/c 2 )( g (r 0 )/g(r GUT )) 1/3 

> 1.4 x 10~ 9 for 10~ 8 Hz</<10 10 Hz. (15) 

Here we have used the normalization in (|l3j ) for G/x/c 2 , Hubble parameter h = 
0.75, and assumed a minimal GUT thermal history. Hence, this lower bound 
is valid up to frequencies / ~ 10 _3 a _1 Hz ~ 10 Hz based on our knowledge 
of the number of relativistic degrees of freedom, g , of the primordial fluid up 
to temperatures T ~ 10 3 GeV. Notice that a measurement of the spectral 
density due to cosmic strings at higher frequencies would sample g at higher 
temperatures. 

We may in turn place a lower bound on the amplitude of the dimensionless 
strain predicted for the gravitational wave emitted by cosmic strings: 

he = 1.3x10- 13 / l y / H gr (/)( T ^-) _1 
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> 3.6 x 10- 23 (-L-) 1 for 10" 6 Hz</<10 8 Hz. (16) 

V1 Hz/ 

The expressions dH-IH) are useful for comparison with the planned sensitivities 
of the forthcoming generation of gravitational wave detectorsEi 

The most promising opportunity to probe for a stochastic gravitational 
wave background due to cosmic strings is through a cross-correlation of the 
observations of the advanced LIGO, VIRGO and LISA interferometers. It is 
estimated that the advanced LIGO detectors will have the sensitivity 


z 3/yi 


= 5.2 x 10 


-25 




1 kHz / 


(17) 


for stochastic waves (equation 125c of Thorne@), sufficient to measure the 


minimum predicted strain (|1 
time. More recent calculation: 


near / ~ 100 Hz in a 1/3-year integration 
confirm that the orientation of the advanced 
LIGO interferometers will be sufficient in order to detect the cosmic string 
gravitational wave background. For the LISA projects, comparison of the 
projected strain sensitivity h c ~ 10 -20 at the frequency / ~ 10~ 3 Hz with (16) 
indicates that the space-based interferometer will be capable of detecting a 
gravitational radiation background produced by a network of cosmic strings. 

Other ground-based interferometric gravitational wave detectors are in de¬ 
velopment or under construction. The GEO600 and TAMA300 detectors, op¬ 
erating near frequencies / ~ 10 3 Hz, may also be capable of measuring a cosmic 
string generated background. 

A network of resonant mass antennae, such as bar and TIGA detectors 
may probe for a stochastic background. Successful detection by these anten¬ 
nae will require improved sensitivity and longer integration time. However, 
cross-correlation between a narrow-band bar and a wide-band interferometric 
detector may improve the opportunities. Estimates of the sensitivity of such a 
system, assuming optimum detector alignments, indicate that 


\Ai„t/ib>2.6 x 10" 19 ^/fi gr (/)( 


( f \ 

3 / 2 / t a b s \ 

V 1 kHz ) 

Vl0 7 sJ 


1/2 


(18) 


is necessary to detect a background S2 gr . Hence, for a 1/3-year observation 
time, the bar and interferometer strain sensitivities at 1 kHz must be better 
than ~ 10~ 23 in order to detect the cosmic string background. 

We stress that the amplitude of the cosmic string gravitational wave back¬ 
ground for frequencies / <) 10 Hz is sensitive to the number of degrees of freedom 
of the cosmological fluid at temperatures T > 10 3 GeV. The amplitude of the 
cosmic string background at LISA-frequencies, near / ~ 10~ 3 Hz, is firm, since 
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Figure 6: Summary of the potential cosmological sources of a stochastic gravitational ra¬ 
diation background, including inflationary models, first order phase transitions and cosmic 
strings, as well as a primordial 0.9K blackbody graviton spectrum (the analogue of the 
blackbody photon radiation). Also plotted are the relevant constraints from the COBE 
measurements, pulsar timings, and the sensitivities of the proposed interferometers. Notice 
that local cosmic strings and strongly first-order phase transitions may produce detectable 
backgrounds in contrast to standard slow-roll inflation models. 


the cosmological fluid near the temperature T ~ 10 MeV is well understood. 
However, at the higher frequencies probed by ground-based detectors, our un¬ 
certainty in the number of degrees of freedom of the cosmological fluid, as 
determined by the correct model of particle physics at that energy scale, may 
reduce the predicted amplitude (|l^) of gravitational radiation. 


7 Conclusion 

We have presented improved calculations of the spectrum of relic gravitational 
waves emitted by cosmic strings. We demonstrated that the effect of a grav¬ 
itational back-reaction on the radiation spectrum of cosmic string loops, for 
which there is an effective mode cut-off n* 1$ 10 2 , may serve to weaken the 
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pulsar timing bound on the cosmic string mass-per-unit-length. Arguing for 
a model of radiation by loops, for which either the spectral index is q > 2 or 
there is an emission mode cut-off n*i$10 2 , we obtain the conservative bound 
G/i/c 2 < 5.4(±1.1) x 10~ 6 due to observations of pulsar timing residuals. We 
believe this bound to be robust, in that the spectrum depends weakly on the 
precise value of the mode cut-off, up to ?r* ~ 10 2 . We have noted the interest¬ 
ing result that the flat, red noise portion of the gravitational wave spectrum 
is sensitive to the thermal history of the cosmological fluid, revealing features 
of the particle physics content at early times. We have also discussed the 
possibility of a low-12 universe and hybrid systems of defects, such as strings 
connected by domain walls or monopoles. Finally, we have pointed out that 
the generation of advanced LIGO, VIRGO and LISA interferometers should 
be capable of detecting the predicted stochastic gravitational wave background 
due to cosmic strings. 

We may place the cosmic string scenario in context with other candL 
date sources for a stochastic background of cosmological gravitational vp.vea23 
in Fig. ^]. These include gravitational waves created during inflationcil and 
through bubble collisions following a first-order phase transitionEa. Viable the¬ 
oretical scenarios within detector sensitivity include strongly first-order phase 
transitions, possibly at the end of inflation, and networks of cosmic strings. At 
this stage, other primordial backgrounds from slow-roll inflation, global topo¬ 
logical defects and the standard electroweak phase transition appear to be out 
of range. The discovery of any of the possible cosmological sources will have 
enormous implications for our understanding of the very early universe and for 
fundamental physics at the highest energies. 


Acknowledgments 

The work of EPSS is supported by PPARC grant GR/K29272. The work of 
RAB is supported by PPARC postdoctoral fellowship grant GR/K94799. The 
work of RRC is supported by the DOE at Penn (DOE-EY-76-C-02-3071). 


References 

1. G. Smoot et al [1990], Ap. J. 360, 685. 

2. T. W. B. Kibble [1976], J. Phys. A9, 1387. Ya. B. Zel’dovich [1980], 
M. N. R. A. S. 192, 663. A. Vilenkin [1981], Phys. Rev. Lett. 46, 1169. 

3. A. Vilenkin & E.P.S. Shellard [1994], Cosmic strings and other topological 
defects (Cambridge University Press). 


19 


4. M.P. McHugh, G. Zalamansky, F. Vernotte & E. Lantz [1996], Phys. Rev. 
D54, 5993. 

5. V.M. Kaspi, J.H. Taylor & M.F. Ryba [1994], Ap. J. 428, 713. 

6. S.E. Thorsett & R.J. Dewey [1996], Phys. Rev. D53, 3468. 

7. R.R. Caldwell, R.A. Battye & E.P.S Shellard [1996], Phys. Rev. D54, 
7146. 

8. R.R. Caldwell and Bruce Allen [1992], Phys. Rev. D45, 3447. 

9. R.A. Battye & E.P.S. Shellard [1994], Nucl. Phys. B423, 260. R.A. 
Battye & E.P.S. Shellard [1995], Phys. Rev. Lett. 75, 4354. R.A. Battye 
& E.P.S. Shellard [1996], Phys. Rev. D53, 1811. 

10. P. Casper [1996], private communication; P. Casper & B. Allen [1995], 
Phys. Rev. D52, 4337. 

11. B. Allen & E.P.S. Shellard [1992], Phys. Rev. D45, 1898. 

12. D. Garfinkle & T. Vachaspati [1987], Phys. Rev. D36, 2229. 

13. R.A. Battye [1995], PhD Thesis (Cambridge University) 

14. T. Vachaspati & A. Vilenkin [1985], Phys. Rev. 31, 3052. 

15. R.A. Battye & E.P.S. Shellard [1997], Work in progress 

16. D. Bennett & F. Bouchet [1990], Phys. Rev. D41, 2408. B. Allen & E. 
P. S. Shellard [1990], Phys. Rev. Lett. 64, 119. 

17. D. Bennett [1986], Phys. Rev. D34, 3592. Erratum: Phys. Rev. D34, 
3932. 

18. C.J.A.P. Martins [1997], astro-ph/9701055. 

19. P.P. Avelino, R.R. Caldwell, and C.J.A.P. Martins [1997], In preparation. 

20. X. Martin & A. Vilenkin [1996], Phys. Rev. Lett. 77, 2879. 

21. X. Martin & A. Vilenkin [1996], gr-qc/9612008. 

22. E.J. Copeland, T.W.B. Kibble & D. Austin [1992], Phys. Rev. D45, 
1000. D. Austin, E.J. Copeland & T.W.B. Kibble [1993], Phys. Rev. 
D48, 5594. D. Austin, E.J. Copeland & T.W.B. Kibble [19951, Phys. 
Rev. D51, 2499. 

23. C.J. Copi D.N. Schramm fe M.S. Turner [1996], FERMILAB-Pub- 
96/122-A, ^stro-ph/960605S . 

24. B. Allen, R. R. Caldwell, E. P. S. Shellard, A. Stebbins & S. Veeraragha- 
van [1996], Phys. Rev. Lett. 77, 3061. 

25. For a survey of planned and present experiments, see Gravitational Wave 
Experiments, ed. E. Coccia, G. Pizzella, and F. Ronga (World Scientific, 
Singapore, 1995). 

26. K. Thorne [1987], in Three Hundred Years of Gravitation, S.W. Hawking 
and W. Israel, eds. (Cambridge University Press). 

27. E.E. Flanagan [1993], Phys. Rev. D48, 2389. 

28. K. Thorne [1994], in Proceedings of the 1994 Snowmass Summer Study: 


20 





Particle and Nuclear Astrophysics and Cosmology in the Next Millenium, 
E.W. Kolb and R.D. Peccei, eds. (World Scientific). 

29. P. Astone, J. Alberto Lobo & B. Schutz [1994], Class. Quant. Grav. 11, 
2093. 

30. R.A. Battye & E.P.S. Shellard [1996], Class. Quant. Grav. 13, A239. 

31. L.P. Grishchuk [1977], Ann. NY. Acad. Sci 302, 439.A.A. Starobinsky 
[1979], JETP Lett. 30, 682. 

32. M.S. Turner & F. Wilczek [1990], Phys. Rev. Lett. 65, 3080. 


21 



